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Within the context of the AdS/CFT correspondence we show that the DC thermoelectric con¬ 
ductivity can be obtained by solving the linearised, time-independent and forced Navier-Stokes 
equations on the black hole horizon for an incompressible and charged fluid. 


INTRODUCTION 


A striking feature of the AdS/CFT correspondence is 
that some fundamental properties of the dual conformal 
field theory (CFT) are captured by the geometry of the 
black hole horizon. The temperature of the CFT is equal 
to the Hawking temperature of the black hole which is 
determined by the surface gravity of the black hole. Sim¬ 
ilarly, the entropy of the CFT is the Bekenstein-Hawking 
entropy which is given by one quarter of the area of the 
black hole event horizon. Although not universal [lj, it is 
also interesting that for a sub-class of holographic black 
holes the shear viscosity, r], is also captured by the area 
of the black hole horizon via r] = s/47r, where s is the 
entropy density ns. 

Here we will argue that, remarkably, the DC thermal 
conductivity and, more generally, the thermoelectric con¬ 
ductivity of the dual field theory, are universally captured 
by physics at the black hole horizon. Specifically, one 
needs to solve linearised, time independent and forced 
Navier-Stokes equations for an incompressible charged 
fluid on the curved horizon. 

The thermal conductivity, a property relevant for all 
dual field theories, determines the heat current, Q ®, 
that is produced after applying a temperature gradient, 
£* = — diT/T , at the level of linear response. If the dual 
field theory has a global 17(1) symmetry, then with the 
additional application of an electric field E l , the heat 
current and electric current, J l , that are produced define 
the thermoelectric conductivities via: 


Q l 


rr" T<\' :t \ / Ej\ 
Ta ij TR ij ) \ Cj ) ' 


( 1 ) 


It is important to emphasise that Q l ,J l are the total 
current fluxes, defined later. In seeking applications of 
the AdS/CFT correspondence to real materials, the DC 
conductivities are important observables to study. Our 
results can be used to determine whether or not the dual 
field theory is a conductor or an insulator and, at a more 
refined level, the temperature dependence of the conduc¬ 
tivities, including the appearance of any scaling laws. 

From a theoretical point of view the DC conductivi¬ 
ties are somewhat subtle to study, however, since they 
are generically infinite unless there is some mechanism 
for momentum to dissipate. A natural framework to 


study momentum dissipation is provided by “holographic 
lattices” Q. Namley, black hole solutions with asymp¬ 
totic behaviour at the AdS boundary that are associated 
with adding sources to the dual CFT which break spa¬ 
tial translations. There has been much interest in these 
black holes since they can realise metal-insulator transi¬ 
tions [H|, @] as well as novel incoherent metals Ml- 

The holographic lattices depend on the holographic ra¬ 
dial direction as well as the spatial directions and hence 
constructing them generically involves solving PDEs. As 
a result most examples that have been studied are one¬ 
dimensional lattices, which break translation invariance 
in just one of the spatial directions. An important excep¬ 
tion is provided by Q-lattices @ (and similar construc¬ 
tions [9|) whose matter content can be used to break 
translational invariance periodically in all spatial dimen¬ 
sions while the metric remains translationally invariant. 

A method for calculating the DC thermoelectric con¬ 
ductivity for Q-lattices and one-dimensional lattices was 
presented in 7 . HO. For these lattices the final re¬ 
sult was expressed explicitly in terms of the black hole 
solution at the horizon. This substantially extended a 
similar result found for the DC electric conductivity at 
zero charge density and with no momentum dissipation 
0 - Here we will show that for generic lattices, breaking 
translations in all spatial directions, one cannot obtain 
such an explicit formula for the DC conductivity. How¬ 
ever, the DC conductivity can always be obtained by 
solving Navier-Stokes equations on the black hole hori¬ 
zon. The earlier results can now be viewed as special 
cases in which the fluid equations can be explicitly solved. 

An early connection between gravity and fluids is the 
membrane paradigm 0. More recently, in the holo¬ 
graphic fluid-gravity correspondence [14{, approximate 
solutions to the gravity equations are obtained by solv¬ 
ing relativistic hydrodynamic equations for the bound¬ 
ary theory via a systematic derivative expansion. The 
Navier-Stokes equations arise after taking a scaling limit 
and hence they too can be captured in a dual gravita¬ 
tional description 00 . These connections are appli¬ 
cable in a hydrodynamic limit. On the other hand in 


17] it was shown how solutions of Navier-Stokes equa¬ 


tions on hypersurfaces in Minkowski space give rise to 
solutions of Einstein’s equations. Here, by contrast, so¬ 
lutions of Navier-Stokes equations on the black hole hori- 
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zon lead to exact transport quantities, given by specific 
two point correlators, in the deformed dual CFT. In ob¬ 
taining these results we do not take a hydrodynamic limit 
of the dual field theory (or any other limit) and our re¬ 
sults apply to arbitrary horizon geometries that arise as 
solutions to the equations of motion. We expect that the 
time-dependent and non-linear generalisation of the fluid 
equations that we obtain will also play a role in studying 
holographic lattices, for example in a suitable hydrody¬ 
namic limit (e.g. see [3 El). We emphasise, though, 
that our results here already show that independently of 
the strength of translation breaking effects and the tem¬ 
perature, a hydrodynamic description of DC transport is 
always possible in terms of a specific fluid living on the 
black hole horizon. 

For simplicity, we focus on holographic lattices of D = 
4 Einstein-Maxwell theory; the main results extend very 
simply to D > 4 as well as the inclusion of other matter 
fields. More details appear in [5oj ]. 


BACKGROUND BLACK HOLES 


We consider the D = 4 bulk action 


m 


S = J d A xV=g (i? + 6—iF 2 ) , (2) 


The unit radius AdS^ vacuum solution is dual to a d = 3 
CFT with a global [/(l)symmetry. We will focus on the 
class of electrically charged, static black holes given by 

ds 2 = — UGdt 2 + — dr 2 + ds 2 (E 2 ), A = a t dt , (3) 

where ds 2 (£ 2 ) = gij(r,x)dx l dx 3 is a metric on a two- 
dimensional manifold, £ 2 , at fixed r. Also, U = U(r), 
while G, F and at are all functions of (r, x l ). 

At the AdSi boundary, as r —> 00 , we have 

[/ —>• r 2 , F 1, G -A G{x), 

9ij (r, x) —t r 2 g i:i (x ), a t (r, x) -A n(x) . (4) 

The spatial dependence of the boundary metric, given by 
G(x), gij{x), corresponds to a source for the stress tensor 
of the dual CFT. Similarly, g(x) is a spatially dependent 
chemical potential for the global 17(1) symmetry. An in¬ 
teresting sub-class of solutions is associated with adding 
spatially periodic sources to a CFT in fiat space. In this 
case the functions are all periodic in the spatial coordi¬ 
nates x l and we can, in effect, take £2 to be a torus. 

The black hole horizon is assumed to be located at r = 
0. By considering the Kruskal coordinate v = t+jjy+• ■ ■ 


we deduce that the near horizon expansions are given by 
U (r) = r ^47 t T + U ^ r + ... ^ , 
at(r,x ) = r (aj 0) (a;)G (0 ) (x) + (x)r + ...^ , 

G(r, x) = G (0) (x) + G (1) (x) r + ... , 

F(r, x) = G (0) (x) + F (1) (x) r + ... , 

9ij (r, x) = gif (x) + gf } {x)r + ... , (5) 

where the factor of G® in the leading term of at{r, x) has 
been added so that electric charge density at the horizon 
is simply \ H = . 


PERTURBING THE BLACK HOLES 


We will consider a perturbation that provides sources 
E, ( for the electric and heat currents, respectively, that 
are linear in t. Specifically, generalising BE El, 
study 


we 


S(ds 2 ) = Sg^dx^dx 1 ' — 2tGUC,idtdx l , 

8A = Sa^dx^ — tEidx 1 + ta t Qdx l , (6) 

with 8g^,Sa^ functions of (r,x l ) while Ei = Et(x), Q = 
Ci(x) are one-forms on £2 and we demand that 

d(Eidx l ) = d((idx l ) = 0 . (7) 


This perturbation solves the time dependence of the 
equations of motion at linear order. 

At the AdSi boundary we demand that the fall-off of 
Sg^ is such that the only sources are parametrised 
by E, C- At the black hole horizon, as r —> 0, regularity 
implies that we must have 

Sgtt = U (5g^ (x) + 0{rfj , Sg rr = ^ (Sglf (x) + 0{rfj , 
Sgij = Sglf (x) + 0(r), 5g tr = 8gif (x) + 0{r), 
6g ti =6glf {x) - &GU ^ + 0(r), 

$gn = ^ (fglf (x) + Q{rfj , 

8a t = 8af (x) + 0(r), Sa r = — (Saf* (a;) + 0(rfj , 

^ a i = + atCi) + 0( r ), ( 8 ) 

with Sglf + Sglf — 2 8gif = 0. Note that the logarithm 
terms combine with the terms linear in time in <[6| . 


Electric and Heat Currents 


We define the bulk electric current density as 

F = V=gF" . (9) 
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At the AdSi boundary we find that J' l \ 00 is the electric 
current density of the dual field theory. The gauge equa¬ 
tions of motion, S7 p FP v = 0, imply 

d.r = o, d r .P = dj (V^gF Ji ) , (10) 

and another equation which won’t play a further role. 

For the heat currents, we want to identify equations 
of motion involving the metric perturbation that have a 
similar structure to the gauge equations of motion. First 
consider a vector k which satisfies 

V At fc' i = 0, V„V (M fc ,,) = ak v , (11) 

for some function a, which would vanish if k is a Killing 
vector. We also write ip = and k IJ F IJU = d v 6 + s„, 

with s a one-form and 6 an globally defined function. In 
the special case that the Lie derivative of F with respect 
to k vanishes we have d[^s^ = 0. We now define the 
two-form G: 

- k^F v ^A a - | {ip - 9) F^ . (12) 

The equations of motion then imply that 

= (a - 6) r + \F v »s p - i {C k FyPA p , (13) 

where is the Lie derivative with respect to k. In our 
case k = dt and at linearised order ip = at + Sat, 6 = 
—at — Sat, s = —Ei dx l + at (i dx l and a / 0 when £ g 0. 

We can now define the following bulk current density 

Q i = V=gG ir . (14) 

From (fl3l) we deduce, in particular, that 

W = 0 , d r Q‘ = -dj (2 V^gGF) . (15) 

By calculating the holographic stress tensor, , we find 

G 3 ' 2 V$t u - MJioo = Qioo - tG 3 / 2 VSt ij Q , (16) 

and we conclude that Q*|oo is the time independent part 
of the heat current density [22^. 

Navier-Stokes on the horizon 

The next step in our analysis is to examine the equa¬ 
tions of motion for the perturbed black holes in the con¬ 
text of a Hamiltonian decomposition with respect to the 
radial direction. As is well known the Hamiltonian is 
simply a sum of constraints. We want to evaluate the 
constraints as an expansion in the radius at the black 
hole horizon. The details of this calculation are techni¬ 
cally involved and will be fully described in j2d]. The 
final results, however, are simple to explain. We find 
that the Gauss law constraint implies that diJ ^ = 0 
where J( 0 -j = J 1 1 H . Furthermore, the t components of 


the momentum constraint, H t . = 0, as well as the Hamil¬ 
tonian constraint, F[ = 0, each separately imply that 
diQ ^ 0 ) = 0 with Q* (0) = Q l \ H - Finally, the i components 
of the momentum constraint, TLj = 0, gives additional 
equations which, when combined with the others, gives 
the linearised Navier-Stokes equations on the black hole 
horizon, presented below. 

To summarise, evaluating the constraints at the hori¬ 
zon leads to a closed system of equations for a subset 
of the perturbation which must be satisfied at the black 
hole horizon. The black hole horizon is as in ([5 ]) -(|5 jl . The 
perturbation at the horizon is given as in (|6jl -(|8 ]> and it 
is illuminating to introduce the following notation: 

Vi = -Sg^\ w = 6a[ 0) , 

A (°) 

P = ~ 47rT - 5(o) v f ln G i0) 6gl ° ] , (17) 

The current densities at the horizon can be written 

J( o) = Phv 1 + erg (djW + Ej) , 

Q(o) = Ts h v i , (18) 

where we define the horizon quantities 

PH = y/g(ojai 0 ' 1 , s H = ^yg{o ), 

4 = V5G)5go)> ^ = ( 19 ) 

The four unknowns in m satisfy the following system 
of four linear partial differential equations: 

vy = 0 , (20) 

V 2 w+t; i V i (ai 0) ) = -V i ^, (21) 

rj H [-2 V*V(j Vj) + Vj p\ = TshCj + Ph (Ej + djw) , 

(22) 

where the covariant derivatives are with respect to the 
metric on the black hole horizon and all indices are 
being raised and lowered with this metric. The first two 
equations are simply diQ y = di Jfa = 0. Note that in 
(E^l) we can also write 2 V I V(jU J -) = V 2 Uj + RjiV ®. 

Remarkably, we have obtained the time-independent, 
linearised Navier-Stokes equations for a forced, incom¬ 
pressible, charged fluid on the curved horizon. Such equa¬ 
tions are also called Stokes equations. The fluid velocity 
is Vi, the effective pressure is p and w is a scalar poten¬ 
tial. The forcing terms are given by the one-forms 47 rT^ 
and E. In (fTHl) and ([22]) . we see that pH , sh, trg and t]h 
can be viewed as coefficients in the constitutive relations 
for the horizon fluid: pn and sh are the charge density 
and the entropy density while erg and t/h are transport 
coefficients associated with the electric conductivity and 
the shear viscosity of the horizon fluid, respectively. We 
stress that 77 #, erg are not , in general, the shear viscos¬ 
ity and DC electric conductivity of the deformed dual 
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CFTf23j]. It is also worth noting that possible thermo¬ 
electric transport coefficients for the black hole horizon 
fluid an , olh and kh are all absent in the expressions 

m 

We now establish a number of interesting properties of 
this set of equations. Firstly, we multiply (E2l) by Vj and 
then integrate over the horizon, leading to 


d 2 x y/g^ 2V {i v j) V {i v j) + (Vw + E) 2 

= j d 2 x (Q{ 0) Ci + J\ 0) E *) . ( 23 ) 


In the case of non-compact horizons we have assumed 
that possible boundary terms vanish. Observe that the 
left hand side is a manifestly positive quantity and this is 
associated with the thermoelectric conductivities being a 
positive semi-definite matrix. 

Second, we consider the issue of uniqueness for (l 20 l) - 
(l 22 l) . If we have two solutions then the difference of the 
functions will satisfy the same equations but with van¬ 
ishing forcing terms, £ = E = 0 . Denoting the difference 
by ( Vi,w,p ), we immediately conclude from ( 1231 ) that 

V ( V)= 0 , V;iz; = 0 . ( 24 ) 


We also have v l diof^ = 0 from (l2lT) and Vp = 0 from 
m- We conclude that the solution space is unique up to 
Killing vectors of the horizon metric, with p , w constant 
and Sg r t fixed by (HB- This results agrees with the in¬ 
tuition that one should be able to boost along the orbits 
of Killing vectors to obtain a solution with momentum. 

Third, we observe that when {E, £) are exact forms, 
{EX) = (de,dz) with e,z globally defined functions on 
£ 2 , we can solve equations m-m by taking w = — e 
and p = 47tTz, plus possible constants, and v 1 = 0. We 
observe that this solution gives zero contribution to the 
current densities (1181) at the horizon. This solution gives 
no contribution to the DC thermoelectric conductivity, 
which we discuss below. Thus, the DC conductivity is 
determined by the harmonic part of E and £. 


The Thermoelectric DC conductivity 

We have shown that the electric and heat currents at 
the horizon, given in (1181) . can be expressed in terms of 
the sources E , £ after solving the Stokes equations (l 20 l) - 
(l 22 l) . To obtain DC conductivities of the field theory we 
need to relate the currents at the black hole horizon to 
the currents at the AdS boundary. In some cases these 
are the same. In general, however, the currents depend 
on the radial coordinate r and one needs to integrate 
V®, OH) - In general, however, we can always define 
total current fluxes which are independent of r and hence 
obtain associated DC conductivities. 


As a concrete example, assume 2J] we have a periodic 
holographic lattice on R 1,2 . That is, the lattice deforma¬ 
tions G{x)Xij{x ) and fi(x) in (j 4 | are periodic functions 
of the spatial coordinates (a; 1 , a; 2 ) ~ (a; 1 + Li,x 2 + L 2 ). 
Defining the total electric current flux densities through 
the x 2 plane or the aq plane, respectively, 


J1 = T 2 . 


J x dx 2 


j 2 = — rdx 1 
Ei J 


( 25 ) 


and defining Q l in a similar way, we can immediately 
deduce from (fTUl) , (fl 5 l) that d r J l = d r Q' = 0, which is 
just Stokes’ theorem in the bulk. These current fluxes at 
the AdS boundary are thus given by their values at the 
black hole horizon which in turn are fixed by E, £ after 
solving (l 20 l) - (l 22 l) to in order to obtain v l and w and then 
using m- This data then gives the DC thermoelectric 
conductivities via AID- 


EXAMPLES 

For the special case that the lattice depends on only 
one of the spatial coordinates we can explicitly solve the 
Stokes equations and recover the results of 11|. We will 
present the details of this calculation in [2o|]. Here we will 
consider the case of a perturbative and periodic lattice 
associated with coherent metals with Drude peaks [25j|. 
Specifically, we consider perturbative solutions about the 
AdS-RN black brane with a flat horizon. If A is the ex¬ 
pansion parameter, at the horizon we assume 

G^ = /(0) + A/(1) H-, gff = gS i] + A h\) -f-, 

^ = 0 -)- A<Z(i) + • • • , (26) 

with /(q), g and a constants, and the remaining functions 
are periodic on the torus £ 2 - For the Ricci tensor we 

have R(g)ij = A + A 2 R^ -\ -and similarly for the 

Christoffel symbols. 

We can solve m-m perturbatively, using the expan¬ 
sion: v = u( 0 ) + j V(i) H-, w = j W( 1 ) + W( 2 ) H- 

and p = A P( 1 ) +P( 2 ) + • • • ■ At leading order we find that 
di i>( 0 ) = 0, = 0 implying that v ^ are constant on 

the torus. At next order we have v = Nfaj with 

n [ 1)3 = -D - 1 [d k {T^y kJ + rWj - d' (n-'djRW)). 

(27) 

Note that the function D -1 / is defined up to a constant 
on a torus. Such constants are fixed at third order in the 
expansion and they do not affect the DC result at leading 
order. 

We next integrate equation (f22l) on the black hole 
horizon, discarding any boundary terms in the non¬ 
compact case, and keep the A 0 part. After defining 
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T,j = (47 tT Q +aEj ), now with Q,Ei constant on the 
torus, we deduce that at leading order in A we have 

v* * {L-yiXj, Jfo^PHV 1 , Q\ 0) ^Tsv\ (28) 

Here pn = g 2 a is the charge density at the horizon, s = 
47 ig 2 is the entropy density, and the constant matrix L 
is given by 

L P = AV 1 J h (^-d i h$d i hW kl + dihft d k N l (1)i 
+ i h (1) djiD-'diRW) + g 2 a {1) d j (□- 1 a i a (1) ) ) , (29) 

where is given by (1271) and 5 l E 

Using (1751) we can write J l « pv l , Q l « Tsv 1 where p = 
(LiL 2 ) _1 / dx 1 dx 2 ^/ zr gF tr \ oa = pn is the total averaged 
charge density. Hence, from m the DC conductivities 
are given by: = (L _1 ) lJ 47rsT, a = a = (L _1 )*- 7 47rp, 

a = (L _1 ) ZJ 47rp 2 / s. 

It is interesting to observe for this general class of 
holographic lattices that at leading order R l] (aT)J h = 
s 2 /p 2 S k , corresponding to a kind of Wiedemann-Franz 
law. Also the thermal conductivity when J = 0, n = 
R — Taa^^a, as well as the electric conductivity when 
<5 = 0, <tq=o = cr — TaR~ x a appear at order A 0 in the 
expansion. These general results complement those using 
other techniques in [2(| . 
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